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Abstract — Synthesizing arrays with low sidelobe and pencil 
beam radiation profile is under investigation for decades. A 
variety of array structures are available, but the simplest and 
useful structure is that of a linear array. Here, two basic 
symmetric Linear Antenna Array structures are assumed. The 
required array structure is assumed to provide low sidelobe 
and pencil beam profile. Departure from a uniformity in 
current and location profile has shown quiet appreciable 
improvement in the radiation pattern. The simulations are 
carried out using Differential Evolution Algorithm employing 
Best of Random mutation strategy (DEBoR). 

Index Terms — Non-uniform Current-Location profile, Low 
Sidelobe-Narrow Beam synthesis, Low Current Tapering, 
Differential Evolution Algorithm,Best of Random Mutation. 

I. Introduction 

From past few decades since the concept of using arrays 
instead of a single element has been evolved, researchers 
have took the challenge to provide various array designs to 
provide a radiation characteristics according to the 
requirements [1-10]. Synthesizing an array depends on several 
matters, like requirements on the radiation pattern, directivity 
pattern etc. Radiation pattern depends on the number and 
the type of elements being used, the physical and electrical 
structure of the array etc. Numerous variations over the 
antenna structures as well as the type of elements are 
available, but for simplicity only one kind of elements are 
used in the whole array structure [11]. Aiming towards a 
radiation pattern with low sidelobe and narrow beam, this 
paper deals with linear array geometries with omnidirectional 
antenna elements without any inter-element phasing. For 
simplicity in formulation of pattern characteristics two 
symmetric structures of linear arrays are assumed, one without 
and the other with a centre element and the whole array is 
assumed to be broadside. The design goal in this paper is to 
find-out a probable profile of locating elements and the current 
distribution over the linear array aperture, that could provide 
low and decreasing sidelobes, narrow main beam, but with 
the least sacrifice in the directional characteristics. In search 
of a good radiation and directional characteristics, non- 
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uniformity in both the current distribution and the inter- 
element distance profile is proposed. The establishment of 
the paper follows Design Equation in section II, A brief 
discussion on Differential Evolution Algorithm in section III, 
Discourse on results in section IV and ends up summing and 
concluding on the results in section V 

II. Design Equations 

Figures 1 and 2 denote two Symmetric Linear Antenna Array 
structures placed along z-axis, without and with a centre 
element respectively. Symmetric linear broadside array 
structure without any centre element, have a generalized Array 
Factor as given by [ 1 1] 



AF w (I,x,6) = 2V /„ cos(Kx n cos 0) 



n=\ 



(1) 



The Array Factor corresponding to a symmetric structure as 
considered in the text with a centre element is given by, 



AF C (I, x, 0) = I + 2 Y /„ cos{Kx„ cos 0) 



n=\ 



(2) 



Where, ]\[ is the number of elements on one side of the array 
axis, 

/ is the excitation amplitude of the n ,h element from the 

centre, and / = |/j , I 2 , ■•■ , I , ■•■ , I N } is the current 
vector defining the current distribution over the aperture. 
Current vector is normalized to its maximum value. I Q is the 
current amplitude for the centre element. 
X„ is the distance of the w th element from the centre, and 



X = \X X , X 2 , ••• , X n , •••, X N j is the location vector. 
Location vector is expressed in terms of the wavelength X ■ 

Ik 

K = ~~ is the wave number, and Q is the azimuth angle. 

A 
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L = 2x N is the length of the array. The inter-element dis- 
tance is assumed never to exceed the limit (/L / 2, X ) . 




Figure 1 . Schematic architecture of a symmetric Linear Array 
Antenna structure of 2N elements placed along z-axis 




Figure 2. Schematic architecture of a symmetric Linear Array 
Antenna structure of 2N+1 elements placed along z-axis 

The directivity of the non-uniform array is given by, 

,2 



D KJ „ = 



Z'< 



"Z"2 T ^ sin{K{z '~ z " )] 



1=1 m=\ 



(3) 



K{z t -z m ) 

Where, N T0T is the total number of the element on the array 

aperture, / ■ and Zj are the current amplitude and location 

of j ,h element on the aperture looking the array from one 
end. 

N T0T = 2N , for a symmetric array without any centre 

element, and 

N T0T = 2N + 1 , for a symmetric linear array with a centre 

element. 
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The Cost Function CF is designed carefully so as to shape 
the optimization as a minimization problem. It is given as 

SLL, 



CF 



SLL C 



+ £|AF(c?J +TR 



+ \BWFNr. - BWFN, 



(4) 



C ""'"!/ 

Where, 

SLL C is the peak sidelobe level (SLL) of the design evolved 

with current generation in dB, 

SLL U is the SLL value in dB of the uniform array with the 

same length as found in the current iteration. 

Z\af(o )l 2 
I \ k A is used to impose null [11] in every 

k 

direction outside the main beam. 

BWFN C denotes the first-null beamwidth (BWFN) of the 

radiation pattern with the parameters generated from current 
iteration, 

BWFN U denotes the BWFN of the radiation pattern of the 

uniform array of same length as that of the non-uniform array 
generated in the current iteration. 

III. The Differential Evolution Algorithm (DEA) 

No Optimization technique has ever been superior on all 
problems. However, according to the survey by . Roscca, G. 
Oliveri and A. Massa [9], since Storn and Price had introduced 
a special type of Evolutionary Algorithm namely, Differential 
Evolution [6], a huge number of researches employing DEA 
have been carried out in the past decade, of which a 
significantly large portion were in the field of electromagnetics 
and antenna synthesis. In [8] a new improved variety of DEA 
is introduced, employing Best of Random selection strategy 
in mutation portion. This is called as Differential Evolution 
with Best of Random algorithm (DEBoR). 

The main idea of DEA is to use vector differences in the 
creation of new probable solutions. DEA aims at evolving a 

population of g trial solutions to achieve a optimal solution 

in K iterations. Each of S individual is identified by one 

chromosome that codes a set of unknown descriptive of the 
problems solution. DEA proceeds in the following steps, 

A. Initialization 

S chromosomes each with JV genes are generated. To 

generate a chromosome, 7 , n — 1,2, ••• N within its lower 

z; in ,n=l,2,---Ar and upper z™, n = 1,2,. ---iV 

bounds a uniform random string r ns G (0, 1) is used in the 
following manner, 
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min . / max min \ 

7 =7 + r \7 — Z I 
Where r varies with n, n = l,2,--- N , 
S, S = 1, 2, • • • S and k, k — 1,2, • • • K m;a . Cost Function 
CF (z) is designed carefully for multi -objective optimization 
process. Maximum number of iterations K , Scaling Factor 
p and Crossover probability CR ls carefully initialized. 
p and CR are random variables in 0.4 < F < 1 and 
0? e (0.3, 0.9]. 
B. Mutation 

3 out of S chromosomes are randomly selected and 
recombined according to the Best of Random (BoR) strategy 
[8] to create a mutant vector v as below, 



.(*) 



F+'lkX-g) 



y is the number of the differential variations, a. , p y and 
Y ' y are the indices of randomly chosen secondary parent 
Z , and two randomly chosen donor vectors ( Zp and 
Z y ), of which, 

CF{zf)< min {CF^) CF^)} Z > is called the 

primary parent. 

C. Crossover 

Binary Crossover is adopted for generating new offspring 

(k) 
solution vector U n s to compete with the corresponding 

primary parent z ■ The process is done as 



Ak) 



Xk) 



Ak) 



if rY:> < CR 
otherwise ' 



Where, r\ n s is a string of uniformly distributed random 

numbers over the range (0, l) . is kept under the same bounds 
as specified for . 

D. Selection 

Based on the fitness is operated on and in order to keep 
better solutions for the next iteration and discard the relatively 
poor one. 



k+\ 



Ak) 



Ak) 



ifCF\u { * ] )<CF(zV 



(k) 



ifCF\z\ K '<CF\u 



(k) 



The steps namely Mutation, Crossover and Selection are 
continued until the maximum iteration is reached. 

IV. Results and Discussions 

This segment establishes the modelled results for different 
broadside symmetric linear antenna array designs optimized 
by DEBoR profiency Eight groups of linear array structures 
are studied. Sets for harmonious structures without core 
elements include 14, 18, 22 and 26 elements, while for centred 
elemented symmetric structures are taken as 15, 19, 23 and 
27. It is found that the best conclusions are received for an 
starting population of 120 chromosomes. Utmost number of 
generations is limited to 400. Scaling Factor p for mutation 
is selected as 0.49. Binary Crossover is chosen with Crossover 
Probability CR as 0.35. 

A. Numerical Data and Resultant Patterns 

Parameters related to the physical, electrical and directional 
properties are Tabulated in Table I, II and III. Inter-element 
distance for a uniform array and Location of elements for a 
non-uniform array is considered as physical parameter, while 
excitation amplitude of each element pair is considered as the 
electrical parameter. Current amplitude for the uniform array 
is considered to be 1. Radiaton parameters are considered as 
the SLL and BWFN. Directivity is commemorated for the non- 
uniform array as radiation parameter. Table I records the SLL 
and BWFN and Peak Directivity of the array sets. The initial 

inter-element spacing is considered as ~ . For symmetric 

linear array sets without and with centre element Table II and 
III record the location of element pairs and the respective 
current amplitudes and compares the SLL and BWFN with 
that of the representing consistent array with same length. 
Peak Directivity as proposed by the optimal non-uniform set 
is marked in each of the Tables . 

Figures 3 and 4 equates the radiation patterns of non- 
uniform array groups to their representing equally long 
uniform array sets for 26 and 27 elements respectively. 

Radiation patterns are diagrammed in linear dB scale in Q - 
plane, since linear arrays as regarded in these papers execute 

omni-directional properties in (j) - plane. 

Table I. SLL, BWFN, and Directivity for Uniform Linear Array Sets 

WITH INTER-ELEMENT SPACING AS /^ / 2 
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Nd. 


Numb ?r o f 
El?m?nts 


SLL 
(dB) 


BWFN 

n 


DtTKtivfe' 


Ll 


14 


-13.30 


16.43 


14 


lb 


15 


-13.32 


::.;. ; 


15 


Ha 


IS 


-13.37 


1276 


IB 


lib 


19 


-13.3E 


12.09 


19 


Ilia 


22 


-13.40 


10.43 


22 


Illb 


23 


-13.40 


9.9S 


23 


IY a 


26 


-13.42 


B.B2 


26 


IYb 


27 


-13.42 


B.53 


-" 
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Table I. Location of Elements, Current Aplitudes, SLL, BWFN, and 

Directivity for Non-Uniform Symmetric Linear Array Sets without 

Centre Element 



Sal 
Nd. 


Location ot 
El*m*nts 


Current 
Amplitudes 


SLL (dB) 


BWFN C) 


Ham 


Opt. 


Unif 


Opt. 


Unif. 


la 


=0.26=0.82 

=1.36=2.05 

=2.87=3.72 

=4.46 


0.90 0.78 

0.96 0.97 

1.00 0.91 

0.75 


-IE. 27 


-13.30 


13.50 


1154 


14.49 


Ha 


=0.2B =0.B2 
=1.42 =2.04 
=2.69 =3.47 
=4.32 =5.16 
=5.90 


0.91 O.EE 
0.S6 1.00 
0.95 0.E5 
0.E4 0.74 
0.69 


-20.63 


-13.37 


11.01 


9. IB 


IE. 63 


:::= 


=0.24 =0.74 
=1.24 =1.75 
=2.25 =2.81 
=3.36 =4.04 
=4.E6 =5.70 
=6.44 


0.7B 0.95 
1.00 0.E9 
0.E3 0.95 
0.92 0.93 
C'.S". C'.~6 
0.73 


-21.13 


-13.40 


10.43 


B.50 


22. B2 


IVa 


=3.28 =0.81 
-1.37 -152 
=2.49 =3.0E 
=3.66 =433 
=4.9E =5.66 
-6.43 -725 


0.97 1.00 
0.91 0.90 
0.90 0.E5 
0.93 0.90 
0.7B 0.70 
0.73 0.75 


-22.02 


-13.42 


B.49 


6.90 


26.97 


=7.99 


0.6E 



Table I. Location of Elements, Current Aplitudes, SLL, BWFN, and 

Directivity for Non-Uniform Symmetric Linear Array Sets with A Centre 

Element 



Set 
No. 


Location of 
El?mints 


Current 
Amplitudes 
h.h.h-ls 


SLL (dB) 


BWFNC) 


D»^-_- 


Opt. 


Unit. 


Opt. 


Unif. 




=0.50 


1.00 0.95 












lb 


=1.01 =153 

=2.17 =254 
=3.77 =4.49 


0.88 0.96 
0.91 0.88 
0.80 0.67 


-19.35 


-13.32 


15.17 


1154 


15.13 




=0.57 


1.00 0.92 














=1.08 =1.63 


0.89 0.88 












lib 


=2.25 =2.85 
=3.54 =437 
=5.18 =5.93 


0.84 0.88 
0.88 0.82 
0.78 0.68 


-20.35 


-13.38 


11.09 


9.16 


1951 




=0.50 


1.00 0.84 














=1 01 =151 


0.88 0.74 












:::i 


=2.02 =2.64 
=3.21 =352 
=4.70 =555 
=6.41 =.11 


0.92 0.93 
0.85 0.91 
0.89 0.77 
0.75 0.67 


-2156 


-13.40 


9 53 


7.71 


23.69 




=0.53 


1.00 0.87 














=1.07 =159 


0.80 0.78 














=2.15 =2.80 


0.90 0.99 












IMj 


=3.56 =427 
=4.86 =555 
=6.34 =7.18 
=8.06 =878 


0.90 0.94 
0.77 0.79 
0.87 0.74 
0.74 0.68 


-20.23 


-13.42 


7.46 


6.29 


27.07 



Figures 5 and 6 portray that an adept selection of position of 
elements and current distribution over the array aperture can 
remarkably mend the resultant radiation pattern, even if 
compared to equally long uniform array. The results are 
tabulated in Tables I, II and III. 

From the Tables I and II and Figure 3 it is clearly viewed 
that, for a 26 element linear array, radiation pattern of both 
the uniform arrays have -13.42 dB whereas, the non-uniform 
array inhibits the sidelobe to -22.02 dB. The Directivity of the 
non-uniform array is 26.97 and that of the initial array is 26. 
The BWFN of the non-uniform array is 8.49°, and those for 
initial and equally long uniform arrays are 8.82° and 6.90°. 



Figure 4 and the Tables I and III show that for 27 element 
array the SLL value of the non-uniform array is reduced to - 
20.23 dB against -13.42 dB of the uniform arrays. Directivity 
of the non-uniform array is 27.07 while that for the initial 
array was 27. BWFN of the optimized set is 7.46°, while those 
for the initial and equally long uniform arrays are 8.50° and 
6.29° respectively. 

R^wkELsr. Pai»ff. sf i B? ^i is lit E yfi-jrje ff Lc 1 ::': :..' r-'\ : f:-7. Array 
with 2fi LtocropLir ElsiTjsr.ct 




■& «a m joa li-a 
Aji^le cf Airival uiDa^ra* 

Figure 3. ©-Plane radiation Patterns of 26 element broadside 

symmetric linear antenna arrays 
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40 BO 1CB 12) 

Af.z!; of Airiv*i iji De Hf «* 

Figure 4. ©-Plane radiation Patterns of 27 element broadside 
symmetric linear antenna arrays 

Ctn=nt E>U tdbt tarn ovar th= tm=ar ipartui= n fctnd v.ithDEA 




^— Location of stftnrjfrr.ts^id - . v^av*t#r.ad-.) ff sjr.cer.tie — ^ 
Figure 5. Optimal Current distribution over the array aperture of 
26 element broadside linear antenna array 
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-s -* -a d a * s 

*£— Location of o to nwrrts^ui-waw length) from centre — ^ 
Figure 6. Optimal Current distribution over the array aperture of 
27 element broadside linear antenna array 

B. Convergence Profile of DEBoR 




ISO 2M ISO 3QQ 

Iteration Cv-rte 



Figure 7. Convergence Profile of Differential Evolution Algorithm 
for optimizing the 26 element broadside symmetric linear array set 

!" 




~ :ef i:iz<r. C-. - : Le 

Figure 8. Convergence Profile of Differential Evolution Algorithm 
for optimizing the 27 element broadside symmetric linear array set 

Figures 7 and 8 registers the convergence profiles of 
Differential Evolution Algorithm as the respective optimization 
processes for 26 and 27 element arrays proceeded. Minimum 

value of the Cost Function CF in eacn iteration is recorded 
for the respective optimization process. These curve display 
that the best up to the current iteration in each process is 
saved and thus results are prohibited from deterioration. The 
optimization processes for 26 and 27 elemental arrays are 
found converging after 337 and 167 iterations. The 
programming has been written in MATLAB language using 
MATLAB 7.5 on core (TM) 2 duo processor, 1 .83 GHz with 2 
GBRAM. 
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Conclusions 

The modelled resultants have demostrated that 
considerable meliorations can be detected with the current- 
location non-uniformity. Less dwindling of currents have 
rendered overall sidelobes in reducing manner and the main- 
beamwidth in between the initial and equally long uniform 
arrays. Thus, it defeats the problems due to unceasing 
sidelobes in a large and the peak directivity is restrained to 
be leastwise to that of the initial array. Hence, Differential 
Evolution with Best of Random mutation scheme has 
executed successfully to work on the linear array synthesis 
problem. 
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